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Abstract 

We study oscillations of accretion torus with azimuthal magnetic field. For 
several lowest-order modes we calculate eigenfrequencies and eigenfunctions and 
calculate corresponding intensity and polarization light curves using advanced ray- 
tracing methods. 



1 Introduction and model description 

In addition to spectroscopy, polarimetry provides us with further information about 
a source of radiation and namely about its geometry. In particular, a time-resolved 
polarimetry may be a useful tool in identifying various types oscillations in non- 
stationary sources. In this note we explore this idea for a system consisting of a 
rotating black hole surrounded by a thick accretion disk (accre tion torus). 

An analytical description of this se t-up was firs t given bv lOkada etafl (|l98fll ) 
in Newtonian gravi t y and recently by Komissarov ( 20061 ) in Kerr geometry (see 



also Montero et al. . 20071 ). We differ from the latter work only by employing a 



polytropic equation of state in the form P = Kp 1+1 / n , but we use the same relation 
between magnetic pressure p m and the enthalpy w, p m = K^C^^w^ with C = 
9ts ~ 9tt94><j> and \x being a parameter (/x = 2 in our work). We assume that the 
angular momentum is constant over the whole volume of the torus. The spatial 
profiles of the fluid density, pressure and the local magnetic field strength measured 
in a local comoving reference frame can be expressed using Lane-Emden functions 
f(r,d) as (Figured]) 

P = P0f n , P = P0f n+1 B = B J^T X f»n ^ 
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where c s is speed of sound and all quantities with subscript zero are evaluated at 
the torus center [ro,7r/2]. At a given point the Lane-Emden function is a solution 
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Figure 1: The view of the torus shape at its cross-section. The plot shows points, where 
there is a maximum of density (p m ), magnetic field induction (B m ) and emissivity (j m ). 
Dashed contours illustrate, where the respective quantities gain one half of their maximal 
values. The solid black line represents the outer boundary of the torus body. 



of the equation 



U - U m + In ^1 + nc s0 f) + 1)( )gJ 2 f] » r J 



(/x-l)(n + l)/3 p [l + nz4/]^ \Co) 

where U = lnu^ is the relativistic effective potential, U m is its value at the inner 
edge of the torus and (3 P = po/PmO is the plasma beta-parameter (/? = 1 in our work). 

2 Polarization from the stationary torus 

Because the fluid is highly magnetized, the dominating emission process in our 
model is synchrotron radiation that produces polarized light. In the local comoving 
reference frame, the direction of polarization is perpendicular to the projection 
of magnetic field onto the polarization plane (i.e. the plane perpendicular to the 
direction of emitted photons). Therefore, in a polarization basis {X, Y}, where 
the Y vector makes an angle ip with that projection, the frequency-integrated local 
Stokes emissivities are given by 

Ji = Jo I — ) I -pr ) ! Jq =pJi cos(20), J v = -pjj sin(2</>), 



PoJ \B 

where Jo is a constant, B± is a normalized projection of the magnetic field onto the 
polarization plane basis and p = (a + l)/(o + 5/3) for the case of the power-law 
distribution of radiating electrons with spectral index a. 

Due to the global orientation of the magnetic field, the synchrotron emission from 
the torus shows certain degree of polarization which depends on the orientation of 
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Figure 2: The dependence of the total luminosity and of the degree and the angle of 
polarization on the location of the observer (right). The left part of the figure shows 
three views of the magnetized torus for inclinations 0°, 45° and 80°. On top of the 
color map of the emission intensity, there are line segments representing the polarization 
properties of radiation at a given place: length of lines means the degree, and orientation 
indicates the angle of polarization. 

the torus with respect to a distant observer. When the observer sees the torus 
from above (along its symmetry axis) the total degree of polarization is zero (or 
negligible) due to axial symmetry. With increasing inclination the net degree of 
polarization increases and is the largest for observers looking edge-on (Figured]). 



3 Polarization changes due to torus oscilla- 
tions 



Accretion tori suffer from various perturbations induced by non-steady accretion 
of mass. Numerical simulations show that perturbations lead to a development of 
variety of oscillating modes in the torus body. In the context of polarization it 
is interesting to examine what effect these internal oscillation have on produced 
synchrotron radiation, because tracing temporal changes in polarization can be a 
possible way to probe the geometry of accretion flows. 

The equations governing linear perturbations of the torus can solved analytically 
in the limit of slender torus, when c s o — > 0. In particular, the torus admits two major 
global modes, radial and vertical epicyclic oscillations, whose poloidal velocity fields 
are nearly uniform on the torus cross-section, and whose eigenfunctions (in terms 
of Eulerian velocity perturbation) can be expressed as 



where m is the azimuthal wavenumber. Frequencies of these oscillations are given 
by co = u r fi + mOo, where CIq and u> r) g is the angular frequency and radial and 
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Figure 3: The dependence of the total luminosity and of the degree and the angle of 
polarization on the time phase of some particular oscillation modes: radial and vertical 
epicyclic modes with wavenumbers m = and m = 1. Lightcurves for three inclination 
are drawn: 10°, 45° and 80°. Horizontal axes show relative phase of each oscillation - 
their absolute periods differ. The degree of polarization for 10° inclination is close to zero, 
hence the corresponding polarization angle curve is not plotted. 
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vertical epicyclic frequencies at the center of the torus, respectively. Recently, the 



radial epicyclic mode has been found also in numerical simulations bv lMontero et al 
(j2007h . 

Using a numerical ray-tracing technique, we calculate the transport of local 
Stokes parameters and construct torus light curves as it would appear to a distant 
observer. We select four distinct oscillation modes, for which we compute time pro- 
files of total intensity, polarization degree and angle. The four modes are radial and 
vertical epicyclic modes with wavenumbers m = and m = l. Resulting lightcurves 
are summarized in Figure [3l Most dramatic changes happen to lightcurves when the 
source is seen under large inclinations (edge-on). That is because Doppler boosting, 
gravitational light bending and other relativistic effects are most pronounced in side 
view. We find that there is a visible difference in temporal evolution of polarization 
lightcurves between radial and vertical modes. There are, however, little differences 
between m = and m = 1 modes (beside some phase shift) of both radial and vertical 
oscillations respectively. Still these modes can be distinguished by their different 
oscillation frequencies. 



4 Conclusions 

In our example, one can discriminate between radial and/or vertical oscillation 
modes easily. On the other hand differences between axisymmetric and non-axisymmetric 
modes are less apparent - the two cases mainly differ in phases but their phase pro- 
files are quite similar. 
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